Unequivocal differentiation of coherent and chaotic hght through interferometric 

photon correlation measurements 
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We present a novel experimental technique that can differentiate unequivocally between chaotic 
light and coherent light with amplitude fluctuations, and thus permits to characterize unambigu- 
ously the output of a laser. This technique consists of measuring the second-order intensity cross- 
correlation at the outputs of an unbalanced Michelson interferometer. It is applied to a chaotic 
light source and to the output of a semiconductor nanolaser whose "standard" intensity correlation 
function above-threshold displays values compatible with a mixture of coherent and chaotic light. 
Our experimental results demonstrate that the output of such lasers is not partially chaotic but is 
indeed a coherent state with amplitude fluctuations. 

PACS numbers: 42.50.Ar, 42.60.Mi, 42.55.Sa 
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The electromagnetic field inside a laser cavity oper- 
ating above the stimulated emission threshold is under- 
stood to be in a quantum mechanical coherent state, 
while below threshold the field is "incoherent" or chaotic, 
consisting of a superposition of waves emanating from 
many independently- and spontaneously-emitting point 
sources [H, [Ij. 

The coherence of the laser output is usually evaluated 
through the second-order intensity autocorrelation func- 
tion of the laser g^^'^lr). Below threshold, g'^^^iO) (i.e. at 
zero time-difference) equals 2, characteristic of a chaotic 
field, while above threshold it takes the value (7^^^(0) = 1, 
characteristic of a quantum coherent state or a classical 
stable wave. When the pump power is varied across the 
threshold in conventional lasers, in which the fraction (/?) 
of spontaneous emission that is coupled to the "useful" 
mode is very small (/? ^ 1), there is an abrupt transition 
between these two values, thus providing a handle for 
the identification of the threshold However, experi- 
ments measuring (7^^^ (0) in high-/? nanolasers have shown 
that g^^^O) undergoes a smooth and gradual transition 
from the value of 2 to 1 0-01 • For pulsed operation, the 
drop is so gradual that a value significantly above 1 per- 
sists even at pumping powers 4 times above threshold 
0. This feature has been interpreted as indicating that 
the output of these nanolasers is a "mixture" of coherent 
(stimulated) and chaotic (spontaneous) light, and does 
not become fully coherent until it approaches the value 
of ff(2)(o) = 1. 

However, according to its definition, .^^^■'(O) measures 
only the variance of the intensity fiuctuations and gives 
no information as to whether the field under test cor- 
responds to a statistical ensemble of randomly-phased 
waves or a coherent quantum state with dynamical am- 
plitude fluctuations. Thus, in conventional low-/? lasers 
the value of g'^^^O) — 1 simply indicates that their output 
has no intensity fluctuations and this arises because their 
gain (i.e. the number of excited dipoles) is clamped. The 
coherent laser field, nevertheless, undergoes phase fluctu- 



ations (corresponding to the Schawlow-Townes phase dif- 
fusion [8|) and these give rise to the flnite spectral width 
of the laser. On the other hand, a coherent laser field 
may also undergo amplitude fluctuations, for example, 
because of noise-excited relaxation oscfllations lol 
In high-/3 lasers, the number of photons at threshold 
(given by 1/ y^) and the number of excited dipoles are 
small, giving rise to "discretization noise" which drives 
relaxation oscillations in the output intensity [l2j and 
causes the value of g'-^' (0) to be signiflcantly higher than 
1 at pumping powers high above threshold 0, [3] ■ Thus 
the value of ^'■^■'(O) > 1 is not necessarily indicative of 
the presence of chaotic spontaneous emission in the laser 
output, but may arise from dynamical fluctuations of the 
coherent field amplitude. 

In this Letter, we present an interferometric photon- 
correlation technique that can unambiguously discrimi- 
nate between coherent and chaotic light. This technique, 
which was initially developed to study spectral diffusion 
in molecules IJ, Il5| , consists of measuring the intensity 
cross-correlation between the two output ports of an un- 
balanced Michelson interferometer in the photon count- 
ing regime, while averaging over many fringes and over 
time. By use of this technique, we demonstrate that the 
output of high-/? nanolasers is indeed a pure coherent 
state (albeit with amplitude fluctuations), and not a sta- 
tistical mixture of coherent and chaotic light. In view 
of this demonstration, we present first some highlights of 
the theory that are necessary for a better understanding 
of the experiments. 

The second-order intensity cross-correlation function 
between the two outputs of a Michelson interferometer 
(see Fig. [IJ can be written as: 

5(2^)(r,<5) = {E%it)E*sit + T)EBit + T)EA{t))t (1) 

where t is the time-difference between the two photon 
detection events, d is the "interferometric delay" corre- 
sponding to the path difference of Sc/n (where c is the 
speed of light in vacuum and n is the refractive index of 
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FIG. 1: Experimental setup. Light is sent into an unbalanced 
Michelson interferometer made of optical fibers. Two single 
photon detectors (SSPD) measure the cross-correlations be- 
tween the two output ports of the interferometer. 



the fiber) between the two interferometer arms, E* /Ei 
is the positive/negative frequency part of the electric 
field operator at output port i = A/B, {...)t denotes a 
quantum mechanical and statistical average, as well as 
an average over time t, while the normalization of g'''^^^ 
through a denominator equal to the product of the two 
output mean intensities is implicit. The output field may 
be written in terms of the input field operators. 
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Em = 
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where a^{t) is the creation operator for the input mode. 
An equivalent pair of equations may be written for 
Ea/Eb and the annihilation operators. Eq. ^ can thus 
be expanded into sixteen terms, among which 10 terms 
average out to zero when the interferometer arm differ- 
ence 6c/n is dithered over a distance of a few wavelengths. 
The six remaining terms are: 

5(2^)(r,,5) = 1 ( a^{S)a\S + T)a{S + T)a{d) 

+ a^O)a\T)a{T)a{0) 

+ {0)a\d + T)a{5 + T)a{0) 

+ a^{S)a\T)a{T)a{S) 

- {S)a'< {T)a{S + T)a{0) 

- a^0)a^6 + r)a{T)a{6)) (3) 

The first two terms correspond each to ^'^^^(t), the 
"standard" second-order autocorrelation function of the 
input field, and arise from the two photons propagating 
like particles along the same interferometer arm, while 
the third and fourth terms correspond respectively to 



g^'^^ (r -|- S) and g^'^^ (r — 5) and arise from the two photons 
propagating like particles along different interferometer 
arms. These four terms are sensitive only to the inten- 
sity fiuctuations of the field. The last two terms, which 
have a negative sign, describe the fiuctuations of the field 
when it undergoes interference by propagating through 
both arms of the interferometer, and are sensitive to both 
amplitude and phase fiuctuations. 

A simple insight into the shape of g^^^^{T,6) may 
be obtained for a strongly unbalanced interferometer 
in which the delay is much greater than the coherence 
time of the incoming radiation {S 3> Tc), so that the 
two peaks corresponding to the displaced autocorrela- 
tion functions ^'■^'(t ± S) are far from the region r w 0. 
Also, for the sake of simplicity, we assume that the field 
presents a Lorentzian spectrum of half- width 1/tc, de- 
termined by both amplitude and phase fiuctuations. We 
denote the correlation time of the amplitude fluctuations 
by r, 
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At r ~ the value of the interference term 
(a^((5)at(0)a((5)a(0)) may be obtained by commuting the 
two creation operators between them, giving to a good 
approximation 



{a\id)a\5)a[5)a{<d)) ^ g'^^\5) « 1 



(4) 



At T > the magnitude of the interference terms 
decreases. For a strongly unbalanced interferometer, 
the correlations dependent on 5 die out, so that for a 
Lorentzian spectrum we have 

{a\5)a\T)a{5 + T)a{Q)) « {a\8)a{5 + T)){a^ {T)a{Q)) 

- |gW(r)p = e-2|^l/- (5) 

For a chaotic field, the second-order correlation function 
obeys the Siegert relation 0| 



g(2)(r)^l + |.g«(r)p 



(6) 



Thus, for a chaotic field the interference terms cancel ex- 
actly the Siegert term of g'^'^\T) and the cross-correlation 
function in the vicinity of zero time-difference is 

5^'aoL(^«0,J) = l (7) 

For a coherent field presenting amplitude fiuctuations, on 
the other hand, the second-order autocorrelation function 
is of the form 
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(2)(t) = 1 + FT[RIN] 



(8) 



where FT[RIN] is the Fourier transform of the Relative 
Intensity Noise (RIN) spectrum of the laser [l^. For a 
Lorentzian RIN spectrum of width 2/Tamp, we have, 



-2|t|/t„ 



(9) 



where a is the integral of the RIN spectrum and cor- 
responds to the variance of the fluctuations divided by 
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the square of the mean intensity. Thus, for a coherent 
field the cross-correlation function in the vicinity of zero 
time-difference is 



■2|T|/r„ 



(10) 



Since < Tamp, the cross-correlation function displays 
a dip at r = which goes down to 0.5(1 + a) and has 
a width of Tc/2. This dip is dug into a broad peak of 
width 2 /Tamp, and this is the signature of a coherent 
field undergoing amplitude fluctuations. We note that 
Eq. ([TU]) differs significantly from what is expected for a 
"mixture" consisting of a fraction x of coherent light and 
(1 — x) chaotic light, which is 



(2X) 
9 mix 



(t«0,<5) = 1 



-2|r|/r, 



(11) 



Our experimental setup is described on Fig. [T] The 
light beam to be tested is introduced in a Michelson inter- 
ferometer set to a path difference oi Sc/n between its two 
arms and equiped with photon counting detectors (su- 
perconducting single photon counters from SCONTEL) 
on both of its output ports. The path difference Sc/ri 
is dithered over several wavelengths, so that the indi- 
vidual fringes at the output of the interferometer are 
time-averaged at the detection process. The output of 
the detectors is thus proportional to the fringe contrast. 
The outputs of the two detectors are fed in a Becker & 
Hickl (B&H) correlation card which plots in real time 
the normalized photon cross-correlation between the two 
output ports g'--^^^{T, S) as function of the time-difference 
(r) between the two photodetection events, with a reso- 
lution of 164 ps. Polarization alignement of the two arms 
is achieved by using Faraday mirrors at the ends of the 
interferometer arms 17|. Blocking one arm of the in- 



terferometer, allows us to measure directly the standard 
second-order autocorrelation function g^'^\T). 

Light from two different sources was studied by means 
of this technique, so as to have a basis of comparison. 
The first source was an unsaturated Erbium Doped Fiber 
Amplifier with no input signal, thus producing spec- 
trally broad Amplified Spontaneous Emission (ASE) with 
chaotic statistics 18]. The ASE was filtered spectrally 
by a a fiber Bragg filter (AOS GmbH) centered at 1559.5 
nm and having a width of 10 pm, in conjunction with 
a spectrally broad Yenista tunable filterjthus producing 
a beam of coherence time Tc ~ 50 ps [19[. The second 
source was a photonic crystal laser operating at room 
temperature and emitting at 1601 nm. The laser cavity 
was formed by eleven missing holes in a line (Lll) in a 
photonic crystal consisting of a triangular lattice of holes 
(diameter 180 nm and period 430 nm) etched in a 260 
nm-thick suspended InP membrane. The gain material 
was a single layer of InAsP quantum dots [20] in the cen- 
tral plane of the membrane. Details on the fabrication of 
this source can be found elsewhere 1211. The nanolaser 



was pumped with a continuous-wave pump at 840 nm, 
and its threshold was determined as the inflection point 
of the logarithmic input-output curve. The laser pre- 
sented a /3 factor of /3 « 1.5 x 10"'^, implying that its 
threshold involves only I/a/S « 25 photons. 
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FIG. 2: Second order autocorrelation function g^^^ir) of a 
nanoscale laser pumped at 1.1 Pth (green dots and curve) 
and a filtered chaotic source (red dots and curve). The points 
are experimental data and the continuous curves are fits to 
Eqs. © and ©. 

The results of the standard g^^^T) measurements on 
both sources, are shown on Fig. [51 The autocorrelation 
function of the chaotic source (red dots) is independent 
of pumping power. At t = it displays a peak whose 
height and width are limited by the experimental resolu- 
tion which is ~ 4 times longer than the coherence time 
(hence g^"^^ (0) = 1.25 rather than 2). The peak can be flt- 
ted by a Gaussian (red continuous curve) with standard 
deviation a « 180 ps, close to the experimental resolu- 
tion. The second-order autocorrelation of the nanolaser, 
on the other hand, depends on pumping power. Its value 
g(2)(o) = 2 below threshold, drops to g'^'^HO) = 1 at 
pumping powers of 2 Pth (with Pth =1.3 mW). The au- 
tocorrelation function of the nanolaser, pumped at 1.1 
Pth, is shown on Fig. [2](green dots). It displays a peak 
of height g^^^(O) = 1.445 which drops to g^'^^oo) = 1 and 
can be fitted by an exponential decay (green continuous 
curve) with time constant Tamp/2 = 1.38 ns, as expected 
from Eq. 

While both light sources display similar second-order 
autocorrelation functions which peak at zero time- 
difference, g'-^HO) > Ij they exhibit a completely dif- 
ferent behavior for their interferometic cross-correlation 
function, g^^"^^ (r, (5), as shown on Fig. [31 For the fil- 
tered ASE source at S = 550 ps, the cross-correlation 
function displays two peaks separated in time by 26, cor- 
responding to two replicas of the autocorrelation function 
displaced to r = ±i5 and reduced by a factor of 4, while 
its value at zero time-difference is g'-'^^\0,6) = 1, which 
is the signature of chaotic light. 
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For the nanoscale laser, on the other hand, the result 
is qualitatively different. The cross-correlation function 
at (5 = 11 ns, in addition to the two replicas of the auto- 
correlation function at r = ±6 (as was the case of chaotic 
light), it displays a strong and narrow dip at zero time- 
difference (r « 0) that drops significantly below 1 and is 
dug into a broad peak with g^^^^r, S) > 1, as expected 
from Eq. (|10|) . This is the signature of a coherent state 
undergoing amplitude fluctuations. If the nanolaser out- 
put were a "mixture" of coherent and chaotic light, the 
dip would lie on a flat background of g^^'^^ (t, (5) = 1 as 
expected from Eq. pT|) . The experimenal curve is well 
fitted by Eq. ([3]) using only the parameters obtained from 
the fit of g''^\T), namely a = 0.445, Tampf^ = 1.38 ns, 
and Tc — 180 ps, the resolution of our setup. 
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FIG. 3: Second order cross-correlation function at the out- 
puts of the interferometer. The points are experimental data 
and the continuous curves are plots of Eq. ^ using the pa- 
rameters extracted from the fits of g^-^^r). (Top) Filtered 
chaotic source. (Bottom) Nanoscale laser pumped at 1.1 Pth- 

The difference between chaotic and coherent light may 
be understood as follows. When a chaotic field is su- 
perposed with itself after a delay much longer than its 
coherence time {6 ^ Tc), as the interfering waves in the 
ensemble have random phases among them, photons may 
come out from either output port of the interferometer, 



giving g(^'^)(0, S) = I. On the other hand, when a coher- 
ent field is superposed with itself after a delay oi 6 Tc, 
even though the value of the phase difference of the two 
parts will be random (implying that the constructive in- 
terference may occur at either one of the two output ports 
of the interferometer), that value will be maintained for 
a time of the order of Tc. Thus, during that time, con- 
structive interference will occur in the same output port 
of the interferometer (photons will come out from the 
same port) , giving rise to an anti-correlation between the 
two ports, that is g'^'^^^{0,S) < 1. In other words, the 
unbalanced interferometer, which superposes two time- 
separated parts of the field, embodies Glauber's view of 
coherence: "In physical optics the term is used to de- 
note a tendency of two values of the field at distantly 
separated points or at greatly separated times to take on 
correlated values... The coherence conditions restrict ran- 
domness of the fields rather than their bandwidth" [l[. 
Thus, measurement of the cross-correlation function at 
zero time-difference gives an unambiguous indicator of 
the coherent or chaotic nature of the field. 

In conclusion, we proposed and demonstrated an ex- 
perimental technique based on the second-order cross- 
correlation of the two outputs of an unbalanced inter- 
formeter, which can unequivocally differentiate a chaotic 
state from a quantum coherent state undergoing dynam- 
ical amplitude fiuctuations. By analysing the output of 
a high-/? nanoscale laser through this technique, we have 
shown that such lasers emit coherent light even when 
their second-order autocorrelation function at zero time- 
difference is greater than unity ^'^•'(O) > 1, resembling 
that of a chaotic state. This value is not due to the 
presence of incoherent spontaneous emission in the laser 
output, but arises from amplitude fluctuations that the 
coherent laser output undergoes during discrete emission 
events, in this system of few emitters and few photons. 
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